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We consider the 3-form field, which has been considered as a candidate for realizing inflation, 
coupled to a scalar field which models the relativistic matter particles produced during the reheating 
epoch. We have investigated the stability conditions for this theory and found that introducing 
£N| ' such a coupling does not lead to any ghosts or Laplacian instabilities. We have also investigated the 

reheating temperature and the production of particles due to parametric resonances. We have found 
that this process is more efficient in this theory compared to the result of the standard-scalar-field 
04 ■ inflationary scenario. 
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I. INTRODUCTION 



The study of cosmology in the last few years has become more and more interesting and fundamental, because 
available data have started giving us a non-trivial picture of the universe at large scales. For example, data led to the 
i-C ■ unexpected result of the acceleration of the universe [3-0|- In the coming years, we may expect new surprises coming 
also in the field of inflation and inflationary non-gaussianities @ . 



In the context of early cosmology, the theoretical predictions have strong connections with fundamental physics. 
However, we still do not have a clear picture of high-energy physics at energies above electroweak scale. Therefore, in 
inflation, many different models have been introduced. People have recently tried to search for alternative models of 
inflation which, on one hand, would still give account for the standard inflationary results, and, on the other hand, 
would leave peculiar and unique-to-the-model imprints in the data. 



Since data are getting more and more precise, some of these models have been already excluded, whereas others are 
still viable. In fact, different models would give different values for each inflationary observable, so that in principle, 
as more data become available, the parameter space in model-space will consistently reduce. 



In the light of high-energy physics and in order to explore the parameter space of models which are not built only 
of fundamental scalar fields, the 3-form inflationary field has been introduced [H, [T3|- Dynamical models for the 3- 
form in cosmological backgrounds have been studied [lll - [l8j . In particular, the background evolution during inflation 
was found to be similar to the one of the standard scalar-field case. However, at the level of perturbations, things 
change considerably [l2l HBj. First, besides the gravitational wave modes — which have standard features — only one 
£> ■ scalar field does propagate. Second, the potential for the 3-form must be chosen such that this scalar mode degree of 
freedom does not become a ghost. Third, the speed of propagation of such a scalar field, is not equal to one, and for 
a class of potentials it may become even negative (leading to a Laplacian instability) . In order to prevent this latter 
possibility to happen, further conditions on the potential must be chosen. Finally, because of this non-trivial speed 
of propagation, these models can — depending on the form of the potential — lead to lar ge v alues of non-gaussianities, 
giving rise to a whole new phenomenological picture for this 3-form inflation theory [1511161 . 

In the standard picture of the early universe, after inflation, the universe passes through the reheating era. It is 
assumed that, at that time, all elementary particles are created and then the universe enters the high temperature 
phase. In order to generate particles during the reheating time, the inflaton field must couple to some matter field 
and, consequently, decay to generate par ticles. This mechanism has been widely studied in the context of modeling 
the matter fields via a scalar field |19l424| (for more recent reviews see e.g. [25J-127J) and f(R) theory [28143 1| . 

In this paper we investigate the role that the 3-form field can play in order to reheat the universe. The model, its 
definition, and the coupling with matter fields are described in Section|TTJ Then, in Section lHTI we study the conditions 
for the theory to be stable against ghosts (i.e. fields with a negative kinetic energy) and Laplacian instabilities 
(i.e. c 2 s < 0) when such a 3-form is coupled with a matter field. Section ITVl is devoted to the study of the background 
approximate solution during the reheating phase. We finally find an estimated expression for the reheating temperature 
and discuss the preheating phase in comparison to other already known and studied models. 
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II. THE LAGRANGIAN 



We will study the Lagrangian for a 3-form field, A^p, |12h15| coupled to a scalar field, </>, given as follows 



S= d 4 x 



rM 2 



pi 



R 
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A „ 1 

(1) 



iv^V^ - \ m 2 2 - ^ E^A aM ^ 



where E a pjg is the Levi-Civita antisymmetric tensor on curved backgrounds, which on Minkowski reduces to e a py$ 
(with eoi23 = 1 = — e 0123 )- Then we also have E a @' yS = e a ^ 1& / ^—g. The last term is a coupling term between the 
3-form and the scalar field (which mimics a relativistic matter field) into which the 3-form decays. This coupling 
term is one of the simplest one we may think of, and it is the lowest dimensional analytical one (in fact here A is 
a dimension- less parameter), which can introduce a decay of the 3-form field into two scalar field particles. Other 
possible coupling terms may arise but we expect their coupling constants to be suppressed by the cutoff scale of the 
theory. By integration by parts, we have that the considered coupling term can be rewritten as 

^ 2 V^, (2) 

where is the vector dual to the 3-form, that is -B M = E^ 01 ^ A a /3 1 /3\. The coupling term can be rewritten as 

_ * E^AafrWut = [-Ai 23 d cj> + A a23 d 1( f> - A O i 3 a 2 + A 012 d 3 <f>} , (3) 

which explicitly shows the coupling between the scalar field and the four independent components of the 3-form. 
We assume a flat Friedmann-Lemaitre-Robertson- Walker (FLRW) manifold 

ds 2 = -dt 2 + a(t) 2 dx 2 , (4) 

and, on this background, we will set the the components of the 3-form A a ^ compatibly with the background sym- 
metries and with the field equations of motion, as in 

Aoij = , Aijk = a 3 eijk X , (5) 
where Cijk is the three-dimensional Levi-Civita symbol (with £123 = 1). The Friedmann equation of motion reads as 

E 1 = 3M^H 2 - Px - P<f> = 0, (6) 
where H = a/a is the Hubble parameter, and 

px = ^X 2 + V + I H 2 X 2 + 3HXX = i Y 2 + V , (7) 

P4> = \<t? + \ ™\<P , (8) 
Y = X + 3HX . (9) 

We also have the second Einstein equation as 

E 2 ee M^(2H + 3H 2 )+ Px + P ^0, (10) 

where px is the 3-form effective pressure defined as follows 

PX = -(^X 2 + V + 3HXX + ^ H 2 X 2 - UVyX 2 ^ . (11) 

P0 = \ft - \ ■ (12) 

The equation of motion for X can be written as 

E x = X + 3{HX + XH) + l2V tV X - A# = . (13) 
Finally, the equation of motion for the scalar matter can be written as 

E ti> ee j> + 3Hj> + ml<P + \<i>Y = . (14) 

The Bianchi identities lead to 

E 1 +3H(E 1 -E 2 ) + YEx+4>Et = 0. (15) 
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III. PERTURBATION THEORY 



In order to study the no-ghost conditions, from the action approach, it is convenient, though not necessary, to 
choose a gauge for which the spatial metric is diagonal (i.e. flat gauge) [36] , so that the metric tensor can be written 

as 

ds 2 = -(1 + 2a)dt 2 + 2d l %pdtdx l +a 2 {l + 2<S>)dx 2 . (16) 
As for the 3-form, we can use a time gauge to fix its scalar perturbations as [l2l . [l5| 

Aoij = ae ijk d k 0(t, x) , A ijk = a 3 e i]k X(t) . (17) 
Finally we will perturb the matter scalar field as <p = <p(t) + Sip. By expanding the action at second order we find 



S 



dta 6 



(d 2 B) 2 (d0) 2 d 2 

^rr~ + w,y^r- + W + uv, y xi> + Y( a + 3$)] 



6V V X 2 



(dip) 2 



i (&M 2 ,H 2 - Y 



■ i d 2 ib 

2M^(Ha - $) - <^J —f 
a 2 - 3M|i$ 2 + - 5(f> 2 



■ 4? 
<9 2 $ 



6M4,iJ$ - 2MS,— 5- + 3(Y 2 + 12V V X 2 )<S> - 
a 2 



mlipi 



Mi 



(a$) 5 



P ^--[12V y X<-Y 



2 ' lMVyyX 4 )® 2 



ipScp — m^cpScp 



1 , 9N r , 1 (dSipf 



(18) 



where we have defined y = A a p 7 A al31 = 6A 2 , on the background. By integrating out the auxiliary fields one finds 
two no-ghost conditions for the remaining two propagating fields ($ and Sip). The independent no-ghost conditions 
for the kinetic matrix can then be written as 



An 



M 2 } a 3 H [6 a 2 (2 HM 2 Y + 3 X 2 H - 2 XY) V y + M 2 x Hk 2 ) 



1 

2 M^H 2 k 2 -3 a 2 (x 2 4> 2 ~ 4 //-' I/, 1 ,, + 4 M^X IIY - (> l/, 2 ,, A'-' II 2 ) \ ' 



'pi 



det A = AuA 22 - A{ 2 = 



3M^a s V y Y 2 



M^H 2 k 2 - 3 a 2 f X 2 <p 2 - 4 H 2 M& + 4 M 2 ,XHY - 6 M 2 ,X 2 H 2 ) y 



pi 



'pi" 



'ph 



Since 



4 H 2 Mp { - 4 M^XHY + 6 M^X 2 H 2 = 
we can change the denominator of the previous two expressions as 

det A = 



3M^a s V y Y 2 



4 -^(i^ + v + Pt ), 



>o, 



M^H 2 k 2 + a 2 V y [AM 2 y (\X 2 + V + P<t> ) - 3X 2 cf> 2 ] 



1 

A 22 = - 



[4M> 2 (i A 2 + V + p^)V y + M^H 2 k 2 } 



2 M^H 2 k 2 + a 2 V y [4M 2 j(iA 2 + V + P4> ) - 3X 2 4> 2 ] 



> 0. 



(19) 
(20) 

(21) 

(22) 
(23) 



These results show that the coupling does not directly contribute to the ghost condition. Condition (1221) . for high fc's, 
is satisfied when 



Vy > , 



(24) 



which corresponds to the no-ghost condition already found in the vacuum case (i.e. in the absence of the (p field) [If 
Furthermore, we now require also 



4M|i(iA 2 + V + P+) - 3X 2 <p 2 = 2M|, (A 2 + 2V) + (2M ] 



\J2 

pi 



3A 



12 



2M 2 x mW 



> 0. 



(25) 
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Sufficient condition for Eq. (|2"5j) to be verified is imposing V > (the other condition found in the vacuum case [If 
and |X/Mpi| < ^/2/3. During reheating (i.e. after inflation ends) this condition is satisfied as, in general, X/Mpi — » 0. 

In order to check whether this condition is violated or not also during inflation, we change the variables to the 
dimensionless variables as 

X 2 Y 2 2 V 2 4> 2 2 m W (0R \ 

w = - ?r, z = - — u = - — v = £ — 7T- 2o 



M P1 ' 6J¥ 2 Af 2, ' 3^^' 6ff 2 M 2 ,' 6^^^ 

Therefore, the no-ghost condition becomes 

(V2 - VSx) 2 + (2 - 3x 2 )u 2 + 2z 2 + 2v 2 > , (27) 

where we have used w 2 ~ 1 which is the requirement for inflation [15| . According to the Friedmann equation ([5]), 
this requirement implies that z 2 <l when the inflation is supposed to be driven only by the 3-form field. Moreover, 
since the contribution of the scalar field <j) to the dynamics of inflationary universe is sub-dominant, the values of 
u 2 and v 2 need to be smaller than z 2 . By taking into account both the first and the second terms, it is clear that 
the no-ghost condition is satisfied when x < \/2/3 and x ^ ^/2/3. However, for x > these terms will take 

a negative value. In order to check that this negative value violates the no-ghost condition or not, one can find the 
minimum of the function f(x), by keeping u constant, as in f(x) = (\/2 — VSx) 2 + (2 — 3x 2 )u 2 . This minimum value 
takes the form / m ; n = — 2u 4 /(l — u 2 ) rj — 2u 4 . Now it is clear that the no-ghost condition is satisfied for all ranges of 
x since / min + 2z 2 + 2i> 2 > 0, as u 2 < z 2 < 1. 



A. Speeds of propagation 

In order to find the two speeds of propagation, we can proceed as follows. We perform the following field redefinition 

* = fcqi, (28) 
A n [6 a 2 (2 M 2 ,H + 3HX 2 - XY) V y + M 2 ,Hk 2 ] 

~A22 qi q2 = H[6a 2 (2 M 2 lH + 3 HX* 2 XY) V, y + Mg.fffc 2 ] k gl + g2 " (29) 
In this case, the second order action, after a few integrations by parts, reduces to 

S = j d t[Q iMi -C^ qj -B ijqi ^ (30) 

where the matrix Qij is diagonal (without approximations, due to the proposed field redefinition) and its diagonal 
elements, for large fc's are only background dependent as they reduce to 

(\n 5 V Y 2 n 3 
Q n ~?±jjgL, Q 13=0l Q 22 = y. (31) 

Along the same lines, by using the equations of motion, we can prove that at order k 2 , the matrix elements CV, reduce 
to 

6a 3 F 2 (V:, + 12^A 2 )fc 2 a* 2 

Wl- ^2 ' { - / 12 — U, 2 2 — -jj - ■ l-J^J 

The antisymmetric matrix By, by using the equations of motion, is of order C(fc _1 ) and can be neglected, as long as 
we consider only the high-fc behavior of the theory. 
Therefore the two speeds of propagation are 

4 = 1+^%^, 4 = 1. (33) 
v ,y 

This result is consistent with the result in vacuum [l5| and with what we would naively expect, due to the form of 
the coupling we have chosen. Since after the end of inflation the ghost-free potentials studied in [l5[ all reduce to 
V ~ ^m 2 X 2 = j^m 2 y, then during reheating there are no Laplacian instabilities as c 2 x — > 1. 
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IV. REHEATING PHASE 



Let us discuss about the behavior of the 3-form field during the reheating phase. For simplicity, we will neglect 
the backreaction effects of the scalar-matter field on the leading term in the background expansion. Because of this 
approximation, we can neglect the interaction term in Eq. (|13p to find a solution of this equation. By differentiating 
Eq. (I13p , one can write this equation in terms of Y as 



Y + 3HY + m 2 Y = 0. 



(34) 



Here, we have used the potential form V = ^rn 2 X 2 for simplicity. This potential form is free from any ghost 
and Laplacian instabilities during inflation [l^|, since the field X is always less than y/2/3Mp\. Moreover, when 
X/Mpi < 1, as already stated above, all potentials which are free from ghosts and Laplacian instabilities at the end 
of inflation approximately take this quadratic form. 

From Eq. (|34p. it follows that Y 2 cx H 2 ^> \H\ when m 2 <C H 2 . so that the three-form field can indeed drive 
inflation. During inflation, H slowly decreases since H < 0, fulfilling in this way the no-ghost condition, up to the 
end of inflation, at scales around m 2 ~ H 2 . After the end of inflation, the reheating phase starts as both the fields X 
and Y begin to oscillate at scales m 2 3> H 2 . The evolution of Y and X during the oscillating phase can be studied 
by changing variable Y — a~ 3 / 2 Y, so that Eq. (13~4"|) becomes 



Y + m 



-H z 



-H)Y = 0. 



Using the approximations m 2 H 2 and m 2 3> \H\, the solution for Eq. (|3"5|) can be written as 

Y = Csin(mi), 

where C is a constant. Thus, the solution of Y can be written as 

Y = Ca~ 3/2 sm(mt). 

By neglecting the contribution from <f> in Eq. © and Eq. (1101) and combining them together, one obtains 



2M, 



pi 



The solution for this equation is 



o(t) 



a n i a 
- + 2 - 

a \a 



c 2 t + 3ci 



= Y 2 = C 2 a~ 3 sin 2 (mt) 



3C 2 [2m 2 t 2 + cos(2mt)] 



16m 2 Af2j 



1/3 



(35) 

(36) 
(37) 

(38) 
(39) 



where c\ and C2 are integration constants. During the oscillating phase, mt 3> 1, so that the scale factor can be 
approximated as 



a(t) 



( 3C 2 t 2 \ 1/3 



By substituting this approximated solution into Eq. (|37l) . one obtains 

Y 



'8 sin(mt) 

— Alp] . 

3 t 

On using the solution for a in Eq. (|40|) . the expression for H becomes 

2m{8c 2 mM| 1 - 3C 2 [sin(2mt) - 2mt]} 



H 



A8m 2 (c 2 t + 3ci)Af 2 i + 9C 2 [2m 2 t 2 + cos(2mi)] 
Therefore, the approximated solution during oscillating phase, mt 3> 1, reads 



H 



M 



(40) 



(41) 



(42) 



(43) 
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Note that this estimated value for H agrees with the numerical result in [ll| since H 2 oc a 3 = t 2 . Substituting H 
into Eq. (| 13[) and on using the relation in Eq. © , the solution for X can be written as 

[sin(mi) — mt cos(mt)] . (44) 

3 {mty 

Substituting X, Y and H into Eq. ([7|), one finds that, at a first approximation, px oc a~ 3 . This implies, as expected, 
that during oscillating phase, the 3- form field behaves as dust. This behavior of the field also agrees with the result 
found in vacuum fla. 



A. Quantum production of particles 

After knowing how the background evolves for the fields in the theory, we now want to find an (approximate) 
expression for the value of the reheating temperature. During the inflationary period, all matter fields are diluted 
away due to the exponentially accelerating expansion of the universe. At the beginning of the oscillating phase, 
all matter fields can then be assumed to start off in their vacuum state. Hence, it is convenient to consider the 
interaction between the classical background field, X, which is driving inflation and the quantum scalar field 4> with 
the Lagrangian in the action ([IJ. The quantum scalar field <fc can be decomposed, due to the symmetries of the FLRW 
manifold, in the Heisenberg picture, as 

*«..>-^/*»M«*~ + <**--). (45) 

where a k and at are the annihilation and creation operators respectively, and k represents the three dimensional wave 
vector. The Fourier mode <j>k{t) obeys the classical equation of motion 

4> k + ZH4> k +(K + ml + XY^ k = 0. (46) 

By comparing this equation of motion to the standard scalar inflaton field ip, with the four-legs coupling interaction 
g 2 ip 2 (p 2 with matter fields, and the three-legs interaction 2g 2 a(f(j) 2 (which arises when the field acquires a vacuum 
expectation value and consequently symmetry gets broken), the modification consists of replacing XY by g 2 (f 2 and 
2g 2 ap respectively, where <p = (2Mpi/v3) sm(m v t) / (m^t) [191 - 124] . On the other hand, as for the f(R) models and 
Starobinsky inflation, defined by f(R) = R + R 2 /(6m 2 ) (32j, the modification in Eq. (|46[) consists of replacing XY by 
£R where R = —4msm(mt)/t [23. l3(j|. It should be pointed out that, considering A as a classical field, implies that 
the coupling term effectively contributes to the (time-dependent) mass term of the modes 4> k as in dm 2 = XY. 



B. Reheating temperature 

In order to find an expression for the reheating temperature in our model, we will use a strategy analogous to the 
one used in /(i?) gravity [29, 30]. First, we introduce a new variable u k = a<fi k and use the conformal time r) = J a _1 dt 
as the time variable. The equation of motion for the modes of the produced particle field can be rewritten as 



d 2 u k 
dr/ 2 



k 2 + a 2 m 2 - -a 2 R + Xa 2 Y 



u k = 0. 



(47) 



where R = 6a~ d a/drj is the Ricci scalar. According to the coupling we have introduced, a heavy particle X at rest 
decays into two particles <fr in which both will have a total energy typically much larger than their rest mass. In other 
words the produced particles will be relativistic, that is (k/a) 2 mi. In this case, we can ignore the second term in 
the parentheses compared to the first one. Since R <x H 2 ~ H , the third term in the parentheses is also, in general, 
negligible when it is compared to the fourth one. In fact, this approximation is valid for large t, since H 2 oc t~ 2 and 
Y oc Mpit -1 . Therefore Eq. (gZJ) will be written as 



d 2 u k 
~d~q 2 ~ 



+ k u k = Uu k , 



(48) 
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where U = —\a 2 Y. For the mode deep inside the Hubble radius, k 2 ^> U, we can choose, as usual, the initial vacuum 
as the state with positive- frequency solution, that is u k = e~ tkr ' /y/2k. The solution of Eq. P%j) can then be written 

as m 

u k = uf + i ^" t/(r/) sin[fc(77 - ?/)K (?/)*/. (49) 

In curved spacetime, the choice of the decomposition of </> into dk and at is not unique. It is possible to use another 
decomposition Af. and A k which can be written in terms of dk and a k as Ak = akdk + P^fik' coefficients otk and 
fik coincide with the coefficients of the Bogoliubov transformation for the ladder operators. This transformation is 
chosen as to diagonalize the Hamiltonian of the field ip at each slicing time rj. Normalization provides the following 
condition for the coefficients ak and 

|a fc | a -|ft| a = l. (50) 

If pk = 0, we get otk = 1 and then we recover the standard creation and annihilation operators. In general, (3 will 
describe the produced particle due to the expansion of the universe and can be written as [34] 



h=-^Cu{rf)e-^drf. (51) 



The energy density of the produced particle in comoving coordinate rj is defined by [33j 

1 f°° 
P V = ^ dkk 2 -k\(3 k \ 2 , 

= A-/ U{ri')dW UMdn dk ke 2lk{r1 ' ~^ . (52) 
8tt 2 Jq Jq Jo 

By using J °° dk ke lkx = — 1 jx 2 and the fact that U both at early and late times, we obtain 

p - = ^L — d7i L —v dv - (53) 

From Y oc a _3 / 2 (i) sin(mt) oc a~ 3 / 2 (rj) sin(J Q v ma dfj), we can estimate U as U = —\a 2 Y oc — Xa 1 ^ 2 (rj) sin(J^ ma dfj) . 
Thus we can write U in terms of conformal time as 

U{rj) = Ca 1/2 {ri) sin( / ' udfj) , (54) 



where C is an overall constant factor and uj = ma. Using the approximation a 1 da/drj <C uj, valid during the 
rapid-oscillations phase, we obtain 

^Ca 1 ' 2 (n)u cos (T udf^j . (55) 
By taking the limit J Q oj dfj 3> 1 and using the relation linifc_ > . 00 sin(kx)/x = nS(x), we obtain 



1 



Pv — ^7T~ / C au} cos / udri drj. (56) 
^ Jo \Jq J 

Note that we have cut out the infinite contribution from our integration. Shifting the phase of the oscillating factor 
by 7r/2 and differentiating the above equation, one obtains 

dp v _ 1 dp v ^ mil 2 ma A \ 2 Y 2 
dt a dr\ 32tt 32tt 

The physical energy density of the matter field relates to the comoving energy density by = p v / 'a 4 . By taking 
into account the effect of the total relativistic matter degrees of freedom, the energy density of the total radiation 
produced in the reheating process takes the following form 

.9* g*m\ 2 /"* 4 2 

fr = ^ = a?' aYdt > (58) 



a 



s 



where t os is a time when the oscillating phase begins. In the regime m(t — t os ) ^> 1, the behavior of the scale factor 
and the averaged expression for Y 2 can be written as 

A A f 2 

a^a (t-t os )^, <y*)^_J^L_. (59) 
Substituting these expressions into Eq. (|551) . we finally obtain 

(60 > 

407T (t - t os ) 

From the evolution of scale factor, one can find the evolution of the Hubble parameter 

ff2 =9ra- (61) 

From Eq. ((BIT)) and Eq. (|6ip. we find that the energy density p r decreases slower than H 2 . Therefore p r becomes, at 
some time, the dominant contribution to the total energy density. We can estimate the duration time of the reheating 
process by using the Friedmann equation ?>M^H 2 ~ p r . This provides (t — t os ) and, consequently the final p r as 

160tt gggAWMgi 
i_t ° 8 - 3 5 *A 2 m' Pr ~ 6400^ 2 ' {b ) 

The energy density of the produced particles is converted to the standard expression for the energy density of a 
radiation gas as in g*ir 2 T^ b /30. Therefore the reheating temperature can be estimated as 

^<A^V /4 M P1 f^V /2 . (63) 



For the f(R) gravity model f(R) = R + R 2 /(6m 2 ), reheating temperature takes instead the form 

By comparing these two results, we find that the reheating temperature of the 3-form field is larger than the one of 
the f(R) gravity model by Mpi/m ~ 10 5 . Note that we estimated the value of 3-form mass, to, by using thepower 
spectrum of the curvature perturbation found in [l^ | and using the observation data from seven-year WMAP (Jj. For 
the simple scalar field model with three-leg interaction, the reheating temperature takes the form 

The result from this last calculation is comparable to the result calculated by using the decay rate obtained from 
particle theory Hl| . 



C. Parametric resonance 



We investigate here the contribution of parametric resonances to the total energy density of produced matter 
particles in the 3-form inflation discussed here. The majority of the energy for the 3-form field, at the end of inflation, 
is stored in the k = mode of the field. When the 3-form field oscillates around the minimum of the potential, its 
energy undergoes coherent oscillations. This dynamics is the same as the one for a standard scalar inflaton field. 
In order to take into account the coherent nature of the 3-form field at the end of inflation, we follow the standard 
procedure by investigating the parametric resonance of the system (22l424j |. We will begin by introducing a new 
variable = a 3 / 2 </>fc. Therefore, Eq. (|4"S"|) can be rewritten as 



ij k +(^ + m 2 +\Y-lH 2 - 3 -H^y k = 0. 



(66) 
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In the limit k 2 ^> H 2 ~ H, one can neglect the last two terms of the above equation. Note that we still neglect the 
backreaction of the 3-form quantum field in the following. Thus, we can still use the solution of Y as in Eq. (|4"Tj) . 
Substituting the solution of Y into Eq. |66|) , we obtain 



y k +[^ + m, + X 7W ^j- ) U = 0. (67) 



The parametric resonance occurs due to the oscillating term. In order to see this behavior, we will introduce a new 
variable, z, defined by m(t — t os ) = 2z — ir/2. Then Eq. (f6"7T) is changed to the Mathieu equation [37| as follow 



d 2 y k 
dz 2 



+ (^U-2gcos(2z))y fc = 0, (68) 



where 



4k 2 4771^ 4V8AMpi 
Ak = —s — tH t- i Q — —p= ■ (69 

a 2 ™ 2 m 2 V3m 2 (t-t os ) V ' 

According to the Mathieu equation, there are instability bands in which the perturbation yt grows exponentially with 
the growth index fj,k — q/2. These instability bands depend on the parameters Ak and q. For a broad resonance, 
these parameters must satisfy the conditions Ak — l 2 and q 3> 1 where I 2 = 1,2,3, .... For narrow resonance, these 
parameter must satisfy the conditions Ak ~ I 2 and q -C 1 where I 2 — 1, 2, 3, .... To guarantee enough efficiency for the 
production of particles, the Mathieu equation's parameters should satisfy the broad-resonance condition. However, 
in general, the parameter q decreases in time. Thus q must take a large enough initial value. 
For the f(R) gravity model, the parameters take the form (3(3| 

4k 2 4m 2 , «£ 



a 2 m 2 m 2 m(t — t os ) 

In order to get a broad resonance, g > 1, the coupling constant £ must take a large value. However, for the 3-form 
model, the coupling constant does not need to be large, since q ~ AMpi/m. 

For a standard scalar inflaton field with four-legs interaction, the parameters take the form [22 

k 2 4m l „ 2g 2 M 2 

A k = ^^ + ^ + 2q, q = / (71) 

am m 3m z (i — t os ) m 

This model naturally gives large values for q, initially. However, it decreases faster than the 3-form model, so that 
the broad resonance tends to disappear more quickly. 

Since the parameter q is initially very large, the field passes through many instability bands. This behavior leads to 
a stochastic change in the growth index fi^ . Therefore, in this situation, we need to analyze the system as a stochastic 
resonance [24j . While the above discussion is based on the assumption that the backreaction is negligible, it will be of 
interest to take into account the effect of backreaction for more realistic model. The exponential growth of the field also 
provides the non-adiabaticity in the change of the frequency co 2 = k 2 / 'a 2 + tti 2 + ^/8/3AAfpi sva.[m(t — t os )]/(t — t os ). It 
should be noticed that, during the parametric resonance regime, the produced particles are far away from equilibrium. 
The study of the thermalization at the end of the parametric resonance regime is also of interest, and we leave it for 
a future project. 

V. CONCLUSIONS 

The inflationary paradigm has joined cosmology with high energy particle physics. In particular, high-energy 
theories may leave their footprint in the data due to peculiar properties of the field which drives inflation. 

In this paper we have investigated the properties of reheating due to the presence of a 3-form field coupled with 
a matter field. We have found that the existence of a coupled matter field does not change the conditions that are 
necessary for ghost-free and Laplacian stabilities. However, to avoid a ghost field to appear, the matter field has 
to be sub-dominant during inflation, but this condition is naturally fulfilled for quite a large variety of dynamics. 
The sub-dominance of a matter field is not required after inflation so that reheating may occur without ghost and 
Laplacian instabilities, leading to a successful production of relativistic particles at a scale 

f 9.9* \ „, r ( m x ' 



T Th < A -^—r Mpi . (72) 
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We have investigated also the production of particles due to parametric resonances. In particular we have found that 
this process can be still modeled by the Mathieu equation. Furthermore, compared to the standard minimally-coupled 
scalar inflaton field scenario, reheating is more efficient, as broad resonances typically survive longer, without the need 
of choosing extremely large values for the coupling constant which models the interaction between the 3-form and 
the (relativistic) ordinary matter fields. Finally, we have shown that indeed, the production of particle at the end of 
inflation, due to the peculiar phenomenology of reheating, gives us a possibility to distinguish experimentally such a 
model from others alternative and viable models of inflation. 
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